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Abstract 

We use a recently obtained resummed quark propagator at finite temperature which takes into 
account both the chromoelectric scale gT and the chromomagnetic scale g^T through the Gribov 
action. The electric scale generates two massive modes whereas the magnetic scale produces a 
new massless spacelike mode in the medium. Moreover, the non-perturbative quark propagator 
is found to contain no discontinuity in contrast to the standard perturbative hard thermal loop 
approach. Using this non-perturbative quark propagator and self-consistent vertices, we compute 
the non-perturbative dilepton rate at vanishing three-momentum at one-loop order. The resulting 
rate has a rich structure at low energies due to the inclusion of the non-perturbative magnetic scale. 
We also calculate the quark number susceptibility, which is related to the conserved quark number 
density fluctuation in the deconfined state. Both the dilepton rate and quark number susceptibility 
are compared with results from lattice quantum chromodynamics and the standard hard thermal 
loop approach. Finally, we discuss how the absence of a discontinuity in the imaginary part of 
the non-perturbative quark propagator makes the results for both dilepton production and quark 
number susceptibility dramatically different from those in perturbative approaches and seemingly 
in conflict with known lattice data. 
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I. INTRODUCTION 


The ongoing ultra-relativistic heavy-ion collision experiments at RHIC and LHC enable 
us to study the quark-gluon plasma (QGP) which is a deconhned state of hadronic matter 
generated at very high temperatures and/or densities. Although the quark-gluon plasma 
may be strongly coupled at low temperatures, at high temperature there is evidence that 
resummed perturbation theory can be used to understand the properties of the QGP. To 
perturbatively study the QGP one needs to have an in-depth understanding of the various 
collective modes. These collective modes can be roughly classihed into three types which are 
associated with different thermal scales, namely the energy (or hard) scale T, electric scale 
gT, and magnetic scale g‘^T, where g is the strong coupling and T is the temperature of the 
system. The majority of studies in the literature have focused on the hard and electric scales, 
since the magnetic scale is related to the difficult non-perturbative physics of confinement. 

Based on the Hard-Thermal-Loop (HTL) resummations [1-3], a reorganization of finite- 
temperature perturbation theory called HTL perturbation theory (HTLpt) was developed 
over a decade ago [4]. HTLpt deals with the intrinsic energy scale T as the hard scale and 
the electric scale gT as the soft scale and has been extensively used to calculate various 
physical quantities associated with the deconfined state of matter. These quantities include: 
the thermodynamic potential and other relevant quantities associated with it [4-23], photon 
production rate [24], dilepton production rate [25, 26], single quark and quark anti-quark 
potentials [27, 28], photon damping rate [29, 30], fermion damping rate [31, 32], gluon 
damping rate [33, 34], plasma instabilities [35-37], jet energy loss [38-43], lepton asymmetry 
during leptogenesis [44, 45], and thermal axion production [46]. 

Although HTLpt seems to work well at a temperature of approximately 2 Tc and above, 
where ~ 160 MeV is the pseudo-critical temperature for the QGP phase transition, 
the time-averaged temperature of the QGP generated at RHIG and LHG energies is quite 
close to Tc- Near the running coupling g is large and the QGP could therefore be 
completely non-perturbative in this vicinity of the phase diagram. In order to make some 
progress at these temperatures, it is necessary to consider the non-perturbative physics 
associated with the QGD magnetic scale in order to assess its role. Unfortunately, the 
magnetic scale is still a challenge for the theoreticians to treat in a systematic manner 
since, although its inclusion eliminates infrared divergences, the physics associated with 
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the magnetic scale remains completely non-perturbative [47]. The fact that the 0{g‘^T) 
correction to the Debye mass receives non-perturbative contributions indicates that the 
background physics is fundamentally non-perturbative [48]. The physics in the magnetic 
sector is described by a dimensionally reduced three-dimensional Yang-Mills theory and the 
non-perturbative nature of the physics in this sector is related with the conhning properties 
of the theory. 

Lattice QCD (LQCD) provides a hrst principles based method that can take into account 
the non-perturbative effects of QCD. Lattice QCD has been used to probe the behavior of 
QCD in the vicinity of T^, where matter undergoes a phase transition from the hadronic 
phase to the deconhned QGP phase. At this point, the QCD thermodynamic functions and 
some other relevant quantities associated with the fluctuations of conserved charges at hnite 
temperature and zero chemical potential have been very reliably computed using LQCD 
(see e.g. [49-56]). In addition, quenched LQCD has also been used to study the structure of 
vector meson correlation functions. Such studies have provided critically needed information 
about the thermal dilepton rate and various transport coefficients at zero momentum [57-60] 
and hnite momentum [61]. 

Calculations in LQCD proceed by evaluating the Euclidean time correlation function 
only for a discrete and hnite set of Euclidean times. To obtain the dilepton rate, one 
needs to perform an analytic continuation of the correlator from discrete Euclidean times to 
reconstruct the vector spectral function in continuous real time. However, this is an ill-posed 
problem. To proceed, the spectral function and hence the dilepton rate in continuous real 
time can be obtained from the correlator in discrete Euclidean times through a probabilistic 
interpretation based on the maximum entropy method (MEM) [62-64], which requires an 
ansatz for the spectral function. Employing a free-held spectral function as an ansatz, the 
spectral function in the quenched approximation of QCD was obtained earlier and found to 
approach zero in the low-energy limit [60]. In the same work, the authors found that the 
lattice dilepton rate approached zero at low invariant masses [60]. In a more recent LQCD 
calculation with larger lattice size, the authors used a Breit-Wigner (BW) form for low- 
energies plus a free-held form for high-energies as their ansatz for the spectral function [57]. 
The low-energy BW form of their ansatz gave a hnite low-energy spectral function and low- 
mass dilepton rate. This indicates that the computation of low-mass dilepton rate in LQCD 
is indeed a difficult task and is also not very clear if there are structures in the low-mass 
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dilepton rate similar to those found in the HTLpt calculation [25]. 

Given the uncertainty associated with lattice computation of dynamical quantities, e.g. 
spectral functions, dilepton rate, and transport coefficients, it is desirable to have an alter¬ 
native approach to include non-perturbative effects that can be handled in a similar way as 
in resummed perturbation theory. A few such approaches are available in the literature: one 
approach is a semi-empirical way to incorporate non-perturbative aspects by introducing a 
gluon condensate^ in combination with the Green functions in momentum space, which has 
been proposed in e.g. Refs. [65-70]. In this approach, the effective n-point functions are 
related by Slavnov-Taylor (ST) identities which contain gluon condensates in the deconfined 
phase as hinted from lattice measurements in pure-glue QGD [71]. The dispersion relations 
with dimension-four gluon condensates in medium exhibits two massive modes [65] (a nor¬ 
mal quark mode and a plasmino mode) similar to HTL quark dispersion relations. This 
feature leads to sharp structures (van Hove singularities, energy gap, etc.) in the dilepton 
production rates [67, 72] at zero momentum, qualitatively similar to the HTLpt rate [25]. 

Using quenched LQGD, Refs. [73, 74] calculated the Landau-gauge quark propagator 
and its corresponding spectral function by employing a two pole ansatz corresponding to 
a normal quark and a plasmino mode following the HTL dispersion relations [25]. In a 
very recent approach [75], a Schwinger-Dyson (SD) equation has been constructed with the 
aforementioned Landau-gauge propagator obtained using quenched LQGD [73, 74] and a 
vertex function related through ST identity. Using this setup the authors computed the 
dilepton rate from the deconfined phase and found that it has the characteristic van-Hove 
singularities but does not have an energy gap. 

In a very recent approach [76] quark propagation in a deconfined medium including both 
electric- and magnetic-mass effects has also been studied by taking into account the non- 
perturbative magnetic screening scale by using the Gribov-Zwanziger (GZ) action [77, 78], 
which regulates the magnetic IR behavior of QGD. Since the gluon propagator with the GZ 
action is IR regulated, this mimics confinement, making the calculations more compatible 
with results of LQGD and functional methods [79]. Interestingly, the resulting HTL-GZ 
quark collective modes consist of two massive modes (a normal quark mode and a plas- 
^ An important aspect of the phase structure of QCD is to understand the effects of different condensates, 
which serve as order parameters of the broken symmetry phase. These condensates are non-perturbative 
in nature and their connection with bulk properties of QCD matter is provided by LQCD. The gluon 
condensate has a potentially substantial impact on the bulk properties, e.g., on the equation of state of 
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QCD matter, compared to the quark condensate. 



mino mode) similar to the standard HTL dispersions along with a new massless spacelike 
excitation which is directly related to the incorporation of the magnetic scale throngh the 
GZ action. This new qnark collective excitation resnlts in a long range correlation in the 
system, which may have important conseqnences for varions physical qnantities relevant for 
the stndy of deconhned QCD matter. In light of this, we wonld like to compnte the dilepton 
prodnction rate and the qnark nnmber snsceptibility (QNS) associated with the conserved 
nnmber flnctnation from the deconhned QGP nsing the non-pertnrbative GZ action. 

This paper is organized as follows. In sec. II we briehy ontline the setnp for qnark 
propagation in a deconhned medinm nsing GZ action. In sec. Ill we caicniate the non- 
pertnrbative dilepton rate and discnss the resnlts. Sec. IV describes the compntation and 
resnlts of non-pertnrbative QNS. In sec. 5 we snmmarize and conclnde. 


II. SETUP 


We know that glnons play an important role in conhnement. In the GZ action [77, 78] the 
issne of conhnement is usnally tackled kinematically with the glnon propagator in covariant 
gange taking the form [77, 78] 

puppl p2 






p2 
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( 1 ) 


where the fonr-momenta P = ^ is the gange parameter, and 7 g is called the Gribov 

parameter. Inclnsion of the term involving 'Jg in the denominator moves the poles of the 
glnon propagator oh the energy axis so that there are no asymptotic glnon modes. Natnrally, 
to maintain the consistency of the theory, these nnphysical poles shonld not be considered 
in the exact correlation fnnctions of gange-invariant qnantities. This snggests that the 
glnons are not physical excitations. In practice, this means that the inclnsion of the Gribov 
parameter resnlts in the ehective conhnement of glnons. 

In QGD, the Gribov ambignity typically resnlts in mnltiple gange-eqnivalent copies and, 
as a resnlt, it renders pertnrbative QGD calcnlations ambiguons. However, the dimensionfnl 
Gribov parameter appearing above can acqnire a well-dehned meaning if the topological 
strnctnre of the SU (3) gange gronp is made to be consistent with the theory. Very recently, 
this has been argned and demonstrated by Kharzeev and Levin [80] by taking into acconnt 
the periodicity of the ^-vacnnm [81] of the theory dne to the compactness of the SU{3) 
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gauge group. The recent work of Kharzeev and Levin indicates that the Gribov term can be 
physically interpreted as the topological susceptibility of pure Yang-Mills theory and that 
conhnement is built into the gluon propagator in Eq. (1), indicating non-propagation of 
color charges at long distances and screening of color charges at long distances in the running 
coupling. This also reconciles the original view Zwanziger had regarding jg being a statistical 
parameter [78]. In practice, jc can be self-consistently determined using a one-loop gap 
equation and at asymptotically high temperatures it takes the following form [76, 82, 83] 


D-1 Nr 


IG = 




( 2 ) 


4y27r 

where D is the dimension of the theory and N^ is the number of colors.^ The one-loop 
running strong coupling, = dvrOs, is 

487r2 


9\T) = 


(33-2iV;)ln(|| 


2 \ ’ 


( 3 ) 


where Nj is the number of quark flavors and Qo is fhe renormalization scale, which is 
usually chosen to be 27rT unless specihed. We £x the scale Aq by requiring that as(l-5 GeV) 
= 0.326, as obtained from lattice measurements [85]. For one-loop running, this procedure 
gives Ao = 176 MeV. 

To study the properties of a hot QGP using (semi-)perturbative methods, the effective 
quark propagator is an essential ingredient. After resummation, the quark propagator can 
be expressed as 


iS-\P) = f-Y.{P), 


( 4 ) 


where S(P) is the quark self energy. One can calculate S using the modified gluon propagator 
(1) in the high-temperature limit to obtain [76] 


S(P) = {igfC, 


oo 

. ^ - K) « -(tgfCr Y. I I 


dfl 

Att 


{K} 


X 


A±(fc,7G) 

AEl 


*70 k ■ 7 *7o - k ■ 7 


iPo k — ^ *Po — k P 4 - — ^ 


± 


( 5 ) 


^ Equation (2) is a one-loop result. In the vacuum, the two-loop result has been determined [84] and the 
Gribov propagator form (1) is unmodified. Only jg itself is modified to take into account the two-loop 
correction. To the best of our knowledge, this would hold also at finite temperature. 
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where ^{k} is a fermionic sum-integral, Sf{K) is the bare quark propagator, and 

± njj +nF(4), 

ei = ± m . ( 6 ) 

where ns and np are Bose-Einstein and Fermi-Dirac distribution functions, respectively. 
The modihed thermal quark mass in presence of the Gribov term can also be obtained as 

OO 

(T'g) = ^ h± (fc, 7 g) . (7) 

_l_ o' di 

0 

Using the modihed quark self energy given in Eq. (5), it is now easy to write down the 
modihed ehective quark propagator in presence of the Gribov term as 

iS'-i(P) = f 4 o 7 o - ■ P, ( 8 ) 


where, keeping the structure typically used within the HTL approximation, Aq and Ag are 
dehned as [76] 


Ao{u,p) = u- j dkkn±{k,'^G) ,p) + ,p)\ , 

2a^C X f 

As{uj,p) =p+ ^ j dkkn±{k,-fG) [Qiiuit,p) + Qi{ut,p)]. 


( 9 ) 


Here the modihed frequencies are dehned as uf = E^{u + k — E\)/k and cuf = E\{u: — k + 
E^)/k. The Legendre functions of the second kind, Qq and Qi, are 

( 10 ) 


n i N _ 1 ^ + P 

Qq[oj,P) = x-ln 


2p oj — p 

Qi{uj,p) = -(1 - uQo{uj,p)). 

p 


( 11 ) 


Using the helicity representation, the modihed ehective fermion propagator can also be 
written as 


•c^DN _ l(7o-7-p) , l(7o + 7-p) 

“ 2 K 2 dI ^ 


( 12 ) 


where D+ are obtained as 


D+{uj,p,'yG) = Ao{(^,p) - As{uJ,p) = u: - p 


2g^Cf 


j dkkn±{k,'yG) 


X [Qoi^l^p) + Ql{d^l,p) + Qoidlt^p) + Qlidjf^p)] , 


T>_(a;,p,7G) = Ao{uj,p) + As{uj,p) = uj + p 


2g^CF 

W 


j dkkn±{k,'yG) 


X [Qoi(^i,p) - QiidJi,p) + Qoiuji,p) - Qi{i^t,p)] ■ (13) 













FIG. 1. Plot of the dispersion relations for different values of 7 ^. In the parenthesis, the first one 
represents a collective excitation mode whereas the second one is the corresponding energy of that 
mode. 

Solving for the zeros of 7 ^) gives the dispersion relations for the collective excita¬ 

tions in the medium. In Fig. 1 we show the resulting dispersion relations for three different 
values of the Gribov parameter 7 ^. In absence of the Gribov term (he, 7 g = 0), there 
are two massive modes corresponding to a normal quark mode with energy a;+ and a 
long wavelength plasmino mode with energy a;_ that quickly approaches free massless 
propagation in the high-momentum limit. These two modes are similar to those found in 
the HTL approximation [25]. With the inclusion of the Gribov term, there is a massless 
mode qg with energy uq, in addition to the two massive modes, and g_ [76]. The extra 
mode go is due to the presence of the magnetic screening scale. This new massless mode is 
lightlike at large momenta.^ In this context, we note that in Ref. [70], such an extra massive 
mode with signihcant spectral width was observed near in presence of dimension-four 
gluon condensates [70] in addition to the usual propagating quark and plasmino modes. 
The existence of this extra mode could affect lattice extractions of the dilepton rate since 
even the most recent LQGD results [73, 74] assumed that there were only two poles (a quark 
mode and a plasmino mode) inspired by the HTL approximation. 

In HTL approximation ( 7 ^ = 0) the propagator contains a discontinuity in complex plane 
^ The slope of the dispersion relation for this extra massless spacelike mode qc exceeds unity in some 
domain of momentum. Thus, the group velocity, doJc/dp, is superluminal for the spacelike mode qq and 
approaches the light cone {duj/dp = 1) from above at high momentum. Since the mode is spacelike, there 
is no causality problem. Instead, this represents anomalous dispersion in the presence of GZ action which 
converts Landau damping into amplification of the spacelike dispersive mode. 






stemming from the logarithmic terms in (13) due to spacelike momentum uP' < ]P. Apart 
from two collective excitations originating from the in-medium dispersion as discussed above, 
there is also a Landau cut contribution in the spectral representation of the propagator due 
to the discontinuity in spacelike momentum. On the other hand, for 7 g 7 ^ 0 the individual 
terms in (13) possess discontinuities at spacelike momentum but canceled out when all terms 
are summed owing the fact that the poles come in complex-conjugate pairs. As a result, 
there is no discontinuity in the complex plane. ^ This results in disappearance of the Landau 
cut contribution in the spectral representation of the propagator in spacelike domain. It 
appears as if the Landau cut contribution in spacelike domain for 7 ^ = 0 is replaced by 
massless spacelike dispersive mode in presence of magnetic scale (70 7 ^ 0). So the spectral 
function corresponding to the propagator for 7 g 7 ^ 0 has only pole contributions. As a 
result, one has 


p±(u;,p) = 2^2(7 g) T ^+) + ^(^ ± ^-) + ^((u ± (Ug)] , (14) 

where has poles at a;+, —cu-, and —oJc and D_ has poles at — 1 ^+, and ujc with a 
prefactor, (o;^ — p"^)/2m‘^{'yc), as the residue. 

At this point we would like to mention that the non-perturbative quark spectral func¬ 
tion obtained using the quark propagator analyzed in the quenched LQCD calculations of 
Refs. [73-75] and utilizing gluon condensates in Refs. [65-67, 70, 72] also forbids a Landau cut 
contribution since the effective quark propagators in these calculations do not contain any 
discontinuities. This stems from the fact that the quark self-energies in Refs. [65-67, 70, 72] 
do not have any imaginary parts whereas in Refs. [73-75] an ansatz of two quasiparticles was 
employed for spectral function based on the LQCD quark propagator analyzed in quenched 
approximation. The spectral function obtained with the Gribov action (14) also possesses 
only pole contributions but no Landau cut. As a result, this approach completely removes 
the quasigluons responsible for the Landau cut that should be present in a high-temperature 
quark-gluon plasma. This is similar to hndings in other nonperturbative approaches [65- 
67, 72-75]. We will return to the consequences of the absence of Landau cut in the results 
and conclusions sections. 

Starting from the Euclidean expression (5), we have nnmerically checked for discontinuities and found 
none. We found some cusp-like structures at complex momenta, but E was found to be C°-continnous 
everywhere in the complex plane. 
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Returning to the problem at hand, the spectral density in (14) at vanishing three mo¬ 
mentum (p = IpI = 0 ) contains three delta function singularities corresponding to the two 
massive modes and one new massless Gribov mode. To proceed, one needs the vertex func¬ 
tions in presence of the Gribov term. These can be determined by explicitly computing 
the hard-loop limit of the vertex function using the Gribov propagator. One can verify, 
after the fact, that the resulting effective quark-gluon vertex function satishes the necessary 
Slavnov-Taylor (ST) identity 


(Fi - F2);.r^(Pi, P2) = S-\P^) - S-\P2) . ( 15 ) 


The temporal and spatial parts of the modihed effective quark-gluon vertex can be written 
as 


r° = og 7 ° + &G 7 • P, 

r = cg 7 * + be P* 7 o + do p* (7 • P), ( 16 ) 


where the coefficients are given by 


Og = 1 
bG = - 


2g^C, 


(2ir) 


W [dkkn±{k,-ia) -i-[iQm + iQm] . 

J OJi — OJ 2 

1 




( 2 vr) 


cg = 1 + 


2g^CF 


V [ dkkn±{k,'yG) -^+ <^< 5 ^ 2 ] , 

^ J uji — u:2 

^ j dkkh±{k,'yG) 


3(cai — oj 2) 


[dQi + 6Qg-6Qt,-6Qt,], 


dG = 


2g^Cf 


Y] / dkkn±{k,-iG) - > 

““ J CUi — 072 


with 


dQnl = QniYl^P) - QniYl^P) for U = 0, 1,2 , 

OJ^I = E'l.{ojm Pk - E^)/k for m = 1 , 2 , 
u:^2 = i^rn — k + F ^)/k ioT m = 1,2 . 

Similarly, the four-point function can be obtained by computing the necessary diagrams in 
the hard-loop limit and it satisfies the following generalized ST identity 

F^r'^"(-Fi, Fi; -F 2 , P 2 ) = r"(Fi - F 2 , -Fp P 2 ) - r"(-Fi - F 2 , Fi; P 2 ) . (17) 
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FIG. 2. The self-energy (left) and tadpole (right) diagrams in one loop order. 


III. ONE-LOOP DILEPTON PRODUCTION WITH THE GRIBOV ACTION 


The dilepton production rate for a dilepton with energy uj and three-momentum q is 
related to the discontinuity of the photon self energy 11^'^ (Q) as [86] 

dR a 11 


;,Discn;(Q). 


( 18 ) 


dud^q 127r^Q^ — 1 27ri 

At one-loop order, the dilepton production rate is related to the two diagrams shown in 
Fig. 2, which can be written as 

d?p 


PO 


(2vr) 


Tr 


S{P) T,{K,Q,-P) SiK) T,{-K,-Q,P) 


Tr 


S{P) TP-P,P;-Q,Q) 


(19) 


where K = P — Q. The second term in (19) is due to the tadpole diagram shown in Fig. 2 
which, in the end, does not contribute since F([ = 0. However, the tadpole diagram is 
essential to satisfy the transversality condition, Q^n'^^(Q) = 0 and thus gauge invariance 
and charge conservation in the system. 

Using the n-point functions computed in sec. II and performing traces, one obtains 

d^p 


n!;(9 = 0) = He^T5^ j 

PO 


X 


(27r)3 
(og + bo) 


+ 


{ac — be)" 


D+ (wi, p, 7 g)^- {(^2 ^ PHg ) D_ (wi, p, 7 g)^+ {(^2 ^ PHg ) 

{cg Pbe + dcY , (cg — beP dcY 


P 


-2cl 


T>+(a;i, p, 7g)F>_ (wa, P, 7g) D- (^^i, P, 1g)D+{uj2, p, 7g) 
1 1 


+ 


D+ (wi, p, 7g)^+ {^2 ,P,1g) D- (^^1, P, 7g)^- {^2, P, 7g) 


( 20 ) 
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The discontinuity can be obtained by Braaten-Pisarski-Yuan (BPY) prescription [25] 
Disc T ^/i(po)/ 2 (go -Po) = 27ri(l - e^‘^) j dui j duj 2 nF{uJi)nF{u 2 ) 

X 5{u — Ui — U 2 ) Pl{uJl)p2{(^2), (21) 


PO 


which, after some work, allows one to determine the dilepton rate at zero three momentum 

D 00 00 

p‘^dp / dui / du2nF{uJi)nF{ui2)d{uj — uii — UJ 2 ) 


dR . lOa^ 1 I 2 

-(g = 0)= 


dud^q 


97r^ 


0 


— 00 —00 


,,2 ,,2\2 

4 ( 1 , ^2 \ 


2puj 


P+(^ 1 ,P)P-(W 2 ,P) 




+ 1 


2puj 

uf + u^- 2p2 - 2m2(7G) . g 


2puj 




( 22 ) 


Using (14) and considering all physically allowed processes by the in-medium dispersion, 
the total contribution can be expressed as 

00 

pp, , lOa^ 1 


dR 


dud^q 


{q = 0 ) = 




p dp X 


0 






/ 


p UJ 


2 2 \ 2 


+ - 2^c) npu,a) ( < 1 


pu J 


K P^ 

+4 - u,, - u,_) (|^) 


/ 


xll- 


2 2 ^ 2 
a;l — cj ' 


2pu 

^^+^^-2p^-2m^(7G) l^ 
2pa; J 


(23) 


Inspecting the arguments of the various energy conserving ^-functions in (23) one can un¬ 
derstand the physical processes originating from the poles of the propagator. The hrst three 
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FIG. 3. Various dilepton processes which originate from the in-medium dispersion with the Gribov 
term. 

terms in (23) correspond to the annihilation processes of q+q+ —)■ 7 *, q-q- —t 7 *, and 
qcqG —t 7 *, respectively. The fourth term corresponds to the annihilation of q+q- —)■ 7 *. On 
the other hand, the hfth term corresponds to a process, —>■ q-j*, where a mode makes 
a transition to a g_ mode along with a virtual photon. These processes involve soft quark 
modes (g+, g_, and qc and their antiparticles) which originate by cutting the self-energy 
diagram in Fig. 2 through the internal lines without a “blob”. The virtual photon, 7 *, in all 
these hve processes decays to lepton pair and can be visualized from the dispersion plot as 
displayed in the Fig. 3. The momentum integration in Eq. (23) can be performed using the 
standard delta function identity 



(24) 


where Xi are the solutions of f{xi) = 0 . 

The contribution of various individual processes to the dilepton production rate in pres¬ 


ence of the Gribov term are displayed in the Fig. 4. Note that in this hgure and in subsequent 


hgures showing the dilepton rate, the vertical axis shows the dimensional late dilepton rate 


dR/d^p = dN/d^xd'^p and the horizontal axis is scaled by the thermal quark mass as to 


make it dimensionless. In Fig. 4 we see that the transition process, g_|_ —)■ q-'y*, begins at 
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FIG. 4. The dilepton production rates corresponding to quasiparticle processes in Fig. 3. 



FIG. 5. Gomparison of dilepton production rates involving various quasiparticle modes with and 
without inclusion of 7 ^. 

the energy a; = 0 and ends up with a van-Hove peak^ where all of the transitions from 
g_i_ branch are directed towards the minimum of the g_ branch. The annihilation process 
® A van-Hove peak [87, 88] appears where the density of states diverges as f'{x)\x=xo = 0 since the density 
of states is inversely proportional to f'{x). 
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FIG. 6. Comparison of various dilepton production rates from the deconfined matter. 

involving the massless spacelike Gribov modes, qcQG 1*, also starts at a; = 0 and falls-off 
very quickly. The annihilation of the two plasmino modes, q-Q- —)■ 7 *, opens up with again 
a van-Hove peak at a; = 2x the minimum energy of the plasmino mode. The contribution of 
this process decreases exponentially. At a; = 2771 ^( 70 ), the annihilation processes involving 
usual quark modes, q+q+ —)■ 7 *, and that of a quark and a plasmino mode, q+q- —>■ 7 *, 
begin. However, the former one {q+q+ —)■ 7 *) grows with the energy and would converge to 
the usual Born rate (leading order perturbative rate) [89] at high mass whereas the later 
one {q+q- —)■ 7 *) initially grows at a very fast rate, but then decreases slowly and hnally 
drops very quickly. The behavior of the latter process can easily be understood from the 
dispersion properties of quark and plasmino mode. Summing up, the total contribution of 
all theses hve processes is displayed in Fig. 5. This is compared with the similar dispersive 
contribution when 7 g = 0 [25], comprising processes q+ — )■ g_ 7 *, q+q+ —t 7 *, q-q- —t 7 * 
and q+q- —)■ 7 *. We note that when 7 ^ = 0, the dilepton rate contains both van-Hove peaks 
and an energy gap [25]. In presence of the Gribov term ( 7 ^ 7 ^ 0), the van-Hove peaks re¬ 
main, but the energy gap disappears due to the annihilation of new massless Gribov modes, 
qcqG -t 7*- 

In Fig. 6 we compare the rates obtained using various approximations: leading-order 
perturbative (Born) rate [89], quenched LQGD rate [57, 60], and with and without the 


16 










Gribov term. The non-perturbative rate with the Gribov term shows important structures 
compared to the Born rate at low energies. But when compared to the total HTLpt rate ® it is 
suppressed in the low mass region due to the absence of Landau cut contribution for 7 g 7 ^ 0. 
It seems as if the higher order Landau cut contribution due to spacelike momenta for 7 g = 0 
is replaced by the soft process involving spacelike Gribov modes in the collective excitations 
for 7 g 7 ^ 0. We also note that the dilepton rate [75] using the spectral function constructed 
with two pole ansatz by analyzing LQGD propagator in quenched approximation [73, 74] 
shows similar structure as found here for 7 g 7 ^ 0. On the other hand, such structure at low 
mass is also expected in the direct computation of dilepton rate from LQGD in quenched 
approximation [57, 60]. However, a smooth variation of the rate was found at low mass. The 
computation of dilepton rate in LQGD involves various intricacies and uncertainties. This 
is because, as noted in sec. I, the spectral function in continuous time is obtained from the 
correlator in hnite set of discrete Euclidean time using a probabilistic MEM method [62- 
64] with a somewhat ad hoc continuous ansatz for the spectral function at low energy and 
also fundamental difficulties in performing the necessary analytic continuation in LQGD. 
Until LQGD overcomes the uncertainties and difficulties in the computation of the vector 
spectral function, one needs to depend, at this juncture, on the prediction of the effective 
approaches for dilepton rate at low mass in particular. We further note that at high-energies 
the rate for both 7 g = 0 and 7 g 7 ^ 0 is higher than the lattice data and Born rate. This 
is a consequence of using the HTL self-energy also at high-energies/momentum where the 
soft-scale approximation breaks down. Nevertheless, the low mass rate obtained here by 
employing the non-perturbative magnetic scale {'Jg 7 ^ 0 ) in addition to the electric scale 
allows for a model-based inclusion of the effect of conhnement and the result has a somewhat 
rich structure at low energy compared to that obtained using only the electric scale ( 7 g = 0 ) 
as well in LQGD. 

We make some general comments concerning the dilepton rate below. If one looks at 
the dispersion plots in Fig. 1 for 7 g = 0, then one hnds that U- falls off exponentially and 
approaches light cone, whereas uj+ does not follow fall off exponentially to light cone, but 
instead behaves as [p -|- m‘^{T)/p] for large p. On the other hand, in the presence of 7 g 7 ^ 0 
both a;_ and ojg approach the light cone very quickly, but again a;+ has a similar asymptotic 
® Since HTL spectral function (i.e, 70 = 0) has both pole and Landau cut contribution, so the HTLpt 
rate [25] contains an additional higher order contribution due to the Landau cut stemming from spacelike 


momenta. 
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behavior as before. This feature of a;+ makes the dilepton rate at large oj in Fig. 6 saturate 
for both 7 g = 0 and 7 g 7 ^ 0 , because the dominant contribution comes from the annihilation 
of two as discussed in Fig. 4. In general, the total dilepton rate in Fig. 6 , behaves as 
~ exp(—(u/T) for 7 g(T) = 0 due to the Landau damping contribution coming from the 
quasigluons in a hot and dense medium. As the Landau cut contribution is missing in the 
7 g(T) 7 ^ 0 case, one hnds a leveling off at low uj. In other words, since the Landau damping 
contribution is absent for 7 g(T) 7 ^ 0 , the rate approaches that of the pole-pole contribution 
for 7 g = 0 as shown in Fig. 5, except in the mass gap region. We further note that the 
LQCD rate [57] matches with Born rate at large u simply because a free spectral function 
has been assumed for large u. On the other hand the LQCD spectral function [57] at low u is 
sensitive to the prior assumptions and, in such a case, the spectral function extracted using 
a MEM [62-64] analyses should be interpreted carefully with a proper error analysis [62]. 
Since the MEM analyses is sensitive to the prior assumption, but is not very sensitive to 
the structure of the spectral function at small ca, the error is expected to be signihcant at 
small bj. The existence of hne structure such as van Hove singularities at small uj cannot be 
excluded based on the LQCD rate [57] at this moment in time. 


IV. ONE-LOOP QUARK NUMBER SUSCEPTIBILITY WITH THE GRIBOV 
ACTION 


We now turn to the computation of the quark number susceptibility (QNS) including the 
Gribov term. The QNS can be interpreted as the response of the conserved quark number 
density, n with inhnitesimal variation in the quark chemical potentials p -|- 5/i. In QCD 
thermodynamics it is dehned as the second order derivative of pressure V with respect to 
quark chemical potential, /i. But again, using the fluctuation-dissipation (FD) theorem, 
the QNS for a given quark flavor can also be dehned from the time-time component of 
the current-current correlator in the vector channel [ 6 , 8 , 90, 91]. The QNS is in general 
expressed as 


-I 


= /? 


d^V 


fi —^0 
00 

f duj 


/i 


ji —^0 
2 


27r 1 — e~^‘^ 


= j d‘^x{Jo{0,x)Jo{0,0)) 
Im noo(a;, 0 ). 


(25) 
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where Jo is the temporal component of the vector current and IIoo is the time-time component 
of the vector correlator or self-energy with external four-momenta Q = (cu,^. The above 
relation in (25) is known as the thermodynamic sum rule [90, 91] where the thermodynamic 
derivative with respect to the external source, p is related to the time-time component of 
static correlation function in the vector channel. 

In order to compute the QNS we need to calculate the imaginary part of the temporal 
component of the two one-loop diagrams given in Fig. 2. The contribution of the self energy 
diagram is 

Tl‘UQ) = N,N,TJ2 f ^'Ti[S{P)T’>(K,Q,-P)S(K)T’>(-K,-Q,P)], (26) 

PO 


where K = P — Q. After performing the traces of the self energy diagram, one obtains 
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where D=p(a;,p, 7 ^) were defined in Eq. (13). We write only those terms of Eq. (27) which 
contain discontinuities 


(og 


ho? 
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Calculating the discontinuity using the BPY prescription given in Eq. (21), one can write 
the imaginary part of Eq. (27) as 


Im n*o = 4WiV/vr(l - 


d^p 

(27r)3 
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Ul 
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(31) 


with 


Cl = Im Il_(a;i,p) = 0, 
C 2 = Ini Zl+(a;i,p) = 0. 


The tadpole part of Fig. 2 can now be written as 

d^p 


n;„« 2 ) = iV/iVoT Y / 
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The four-point function Too at zero three-momentum can be obtained using Eq. (17) giving 

(34) 
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Proceeding in a similar way as the self-energy diagram, the contribution from the tadpole 
diagram is 


Im n*o = -4iVciV/7r(l - 


d^p 

(27r)3 
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Cip^{uJ2,p) + C2P^{UJ2,P) 
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dui J du 2 dipj — uji— 0J2) 
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(35) 


The total imaginary contribution of the temporal part shown in Fig. 2 can now be written 
as 


Im IIoo = Im Hgo -h Im Uqq 
= 4:NcNf7i{l - e^‘^) 
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duji / du2 6{u — ui — U2)nF{uJi)nF{uj2) 
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00 
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(36) 


It is clear from (31) and (35) that the tadpole contribution in (35) exactly cancels with the 
second term of (31) even if Ci and C 2 are hnite, e.g., for the HTL case ( 7 ^ = 0) [6, 8]. At 
hnite 7 ^, the form of the sum of self-energy and tadpole diagrams remains the same, even 
though the individual contribution are modihed. 

Putting this in the expression for the QNS in Eq. (25), we obtain 


X,(T) = 


d^p 

(27r)3 


du / dui / di02 5 {oo — ooi — oo2)nF{ooi)nF{oo2) 
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1 - 1 P+(^^1,P)P-(^^2,P) 


.2 \ 2 


= 4iV,iVy/3 




= x7(n 


(37) 


where we represent the total Xq{T) as Xq^{T) since there is only the pole-pole contribution 
for 7 g 7 ^ 0. However for 7 g = 0 there will be pole-cut (xP^(T)) and cut-cut (x“(T)) 
contribution in addition to pole-pole contribution because the spectral function contains 
pole part -|- Landau cut contribution of the quark propagator. 

In Fig. 7 we have presented the different contribution of QNS scaled with the corre¬ 
sponding free values with and without the Gribov term. We, at hrst, note that the running 
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FIG. 7. QNS scaled with free values are compared with and without the inclusion of jg- In each 
case a band appears due to the choice of the two renormalization scales as 27rT and dvrT. The 
various symbols correspond to LQCD data from various groups labeled as WB [49], BNL-BI(B) 
and BNL-BI(u) [51, 52], and TIFR [56]. 

coupling in (3) is a smooth function of T around and below Tc. We have extended to low 
temperatures as an extrapolation of our high-temperature result even though our treatment 
is strictly not valid below T^. Now from the hrst panel of Fig. 7, the pole-pole contribution 
to the QNS with the Gribov action is increased at low T, compared to that in absence of the 
Gribov term. This improvement at low T is solely due to the presence of the non-perturbative 
Gribov mode in the collective excitations. However, at high T both contributions become 
almost same as the Gribov mode disappears. There are no pole-cut (pc) or cut-cut (cc) 
contribution for 7 g(F) ^ 0, compared to that for = 0. The pc and cc contributions in 
absence of magnetic scale are displayed in second and third panels. As a result, we hnd that 
the QNS in presence of magnetic scale contains only the pp-contribution due to collective 
excitations originating from the in-medium dispersion whereas, in absence of magnetic scale. 
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the QNS is enhanced due to additional higher order Landau cut (i.e., pole-cut -|- cut-cut) 
contribution as shown in the fourth panel. When compared with LQCD data from various 
groups [49, 51, 52, 56], the QNS in presence of magnetic scale lies around (10 — 15)% below 
the LQCD results whereas that in absence of magnetic scale is very close to LQCD data. 
This is expected due to the additional higher-order Landau cut contribution in absence of 
magnetic scale as discussed earlier. This also suggests that it is necessary to include higher 
loop orders in QNS in presence of the magnetic scale, which is beyond scope of this paper. 
However, we hope to carry out this non-trivial task in near future. 


V. CONCLUSIONS AND OUTLOOK 

In this paper we considered the effect of inclusion of magnetic screening in the context of 
the Gribov-Zwanziger picture of conhnement. In covariant gauge, this was accomplished by 
adding a masslike parameter, the Gribov parameter, to the bare gluon propagator resulting 
in the non-propagation of gluonic modes. Following Ref. [76] we obtained the resummed 
quark propagator taking into account the Gribov parameter. A new key feature of the re¬ 
sulting resummed quark propagator is that it contains no discontinuities. In the standard 
perturbative hard-thermal loop approach there are discontinuities at spacelike momentum 
associated with Landau damping which seem to be absent in the GZ-HTL approach. Using 
the resulting quark propagator, we evaluated the spectral function, hnding that it only con¬ 
tains poles for 7 g % 0. We then used these results to compute (1) the dilepton production 
rate at vanishing three-momentum and (2) the quark number susceptibility. For the dilepton 
production rate, we found that, due to the absence of Landau damping for 7 g % 0, the rate 
contains sharp structures, e.g. Van Hove singularities, which don’t seem to be present in 
the lattice data. That being said, since the lattice calculations used a perturbative ansatz 
for the spectral function when performing their MEM analysis [62] of the spectral function, 
it is unclear how changing the underlying prior assumptions about the spectral function 
would affect the hnal lattice results. Moreover, the error analsis for spectral function with 
MEM prescription [63] has to be done carefully than it was done in LQGD calculation [57]. 
Since the result is sensitive to the prior assumptions, the error seems to become large and 
as a result no conclusion can be drawn for hne structures at low mass dileptons from the 
LQGD result. For the quark number susceptibilities, we found that, again due to the ab- 
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sence of Landau damping for 7 ^ 7 ^ 0, the results do not agree well with available lattice 
data. This can be contrasted with a standard HTLpt calculation, which seems to describe 
the lattice data quite well with no free parameters. It is possible that higher-order loop cal¬ 
culations could improve the agreement between the Gribov-scenario results and the lattice 
data; however, the success of HTLpt compared to lattice data as well as nonperturbative 
model calculations suggests that at T > 200 MeV the electric sector alone provides an ac¬ 
curate description of QGP thermodynamics. Nevertheless, the present HTLpt results poses 
a serious challenge to the Gribov scenario for only inclusion of magnetic mass effects in the 
QGP. The absence of quasigluons responsible for the Landau cut makes the results for both 
dilepton production and quark number susceptibility dramatically different from those in 
perturbative approaches. We conclude that the results with present GZ action is in conflict 
with those in perturbative approaches due to the absence of the Landau cut contribution in 
the non-perturbative quark propagator. 
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